Nonlocal observables and light cone-averaging in relativistic thermodynamics 



Jorn DunkcEl 

Rudolf Peierls Centre for Theoretical Physics, University of Oxford, 
1 Keble Road, Oxford 0X1 3NP, United Kingdom 



o 
o 

(N 



(N 
o 

S 



I 

o 



(N 
O 

o 



X 
S3 



Peter Hanggi 

Institut fiir Physik, Universitdt Augsburg, Universitdtsstrafie 1, D-86135 Augsburg, Germany 



Argelander-Institut fiir Astronomic, 



Stefan Hilbert 

Universitdt Bonn, Auf dem Hiigel 71, D-53121 Bonn, Germany 
(Dated: March 25, 2009) 



The unification of relativity and thermodynamics has been a subject of considerable debate over 
the last 100 years. The reasons for this are twofold; (i) Thermodynamic variables are nonlocal 
quantities and, thus, single out a preferred class of hyperplanes in spacetime. (ii) There exist dif- 
ferent, seemingly equally plausible ways of defining heat and work in relativistic systems. These 
ambiguities led, for example, to various proposals for the Lorentz transformation law of temperature. 
Traditional 'isochronous' formulations of relativistic thermodynamics are neither theoretically sat- 
isfactory nor experimentally feasible. Here, we demonstrate how these deficiencies can be resolved 
by defining thermodynamic quantities with respect to the backward-lightcone of an observation 
event. This approach yields novel, testable predictions and allows for a straightforward-extension 
of thermodynamics to General Relativity. Our theoretical considerations are illustrated through 
three-dimensional relativistic many-body simulations. 



PACS numbers: 



Introduction.- Thermodynamics, in the traditional 
sense, aims at describing the state of a macroscopic sys- 
tem by means of a few characteristic parameters {Zi\ [l|, 
0, B Typical candidates for thermodynamic state 
variables {Zi\ are either conserved (extensive) quanti- 
ties, e.g. the particle number N and internal energy U , 
or external control parameters that quantify the break- 
ing of symmetries [l[ . Examples of the latter include the 
volume y of a confining vessel, indicating the violation 
of translational invariance, or external magnetic fields, 
which may break the spatial isotropy. Each extensive 
state variable is accompanied by an intensive quantity 
Zi ~ d$ jZi , derived from a suitably defined entropy func- 
tion(al) §({Z.i}). Representing an abstract mathematical 
theory of differential forms [J], thermodynamic concepts 
have been successfully applied to vastly different areas, 
ranging from microscopic many-particle systems 0, 0, @ , 
where § is usually interpreted as an information mea- 
sure (canonical ensemble) or integrated phase space vol- 
ume (microcanonical ensemble), to exotic objects such 
as black holes where § is related to the black hole's 
surface area. 

As a coarse-grained macroscopic theory, thermody- 
namics is inherently nonlocal in that it only considers 
certain global, or averaged, properties of a physical sys- 
tem 0, 3 . This is rather unproblematic within non- 
relativistic Newtonian physics, where statements such as 
"the total energy of a system at time i" are unambigu- 
ously defined for arbitrary observers. By contrast - owing 



to the absence of a universal time parameter - the nonlo- 
cal character of thermodynamics has stirred considerable 
confusion 0, i, i, [p, |lU E^M Q [11 in Einstein's 
theory of relativity [1^, [3 ■ 

To illustrate the conceptual difficulties in relativis- 
tic thermodynamics, consider a confined gas described 
by a particle current density j^{t,x) and an energy- 
momentum tensor density 9^'^{t,x). If the gas is sta- 
tionary in some inertial frame E, then j'^ is conserved, 
i.e., = 0, but the divergence of 9^'^ does not identi- 
cally vanish (due to the pressure arising from the spatial 
confinement, see example below): 



1,2,3. 



(1) 



This means that space-like surface integrals over are 
independent of the underlying three-dimensional hyper- 
surface fi in (l-|-3)-dimensional Minkowski spacetime M4, 
whereas those over B^" do depend on f^. The latter fact 
is problematic since thermodynamic state variables like 
energy or momentum U = (U"'^,U^,'U^) are usually 
defined as surface integrals over the energy- momentum 
tensor (see App.[^and |B]for notation) [Iq : 



M,i^-0,1,2,3. (2) 
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Hence, the first task in relativistic thermodynamics is 
to identify those hypersurfaces that are suitable for 
defining state variables. Subsequently, one still needs to 
settle for appropriate definitions of entropy, heat, etc. 

We shall begin by reviewing how these problems 
are tackled in the most popular, competing versions 
of relativistic thermodynamics, originally proposed by 
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Planck 01 and Einstein [1], Ott and Van Kampen [1], 
respectively. A careful analysis elucidates that the tra- 
ditional approaches are neither conceptually satisfactory 
nor experimentally feasible. The deficiencies can be cured 
by defining thermodynamic quantities in terms of light- 
cone averages. To clarify these as pec ts, we consider a 
weakly interacting relativistic gas [201. Notwithstand- 
ing, the main conclusions apply to any confined system 
that can be described by tensor densities , .... In 

the second part, we shall discuss observable consequences 
such as the apparent drift of distant objects that are, in 
fact, at rest relative to the observer. This surprising ef- 
fect - which should be accounted for when estimating the 
velocities of very hot astrophysical objects from photo- 
graphic data - will be illustrated by relativistic many- 
particle simulations. 

Model (Jiittner gas).- We consider an enclosed, dilute 
gas consisting of N relativistic particles [rest mass m; 
velocity v; momentum p = mt)(l — speed of light 
c = 1]. Let us assume the gas is stationary in the ('lab'-)- 
frame S, and can be described by a S-time- independent, 
normalized one-particle phase space probability density 
function (PDF) 

/(t,a;,p) = (p(a3,p) = £>(a3) </)j(p), (3a) 

with Jiittner momentum distribution [20l . [2ll | 

h{p) = Z"' exp(-/3/), /3 > 0. (3b) 

Z = Ai:m^ K2{l3m) / {f3m) is the normalization constant, 
Kn{z) the nth modified Besscl function of the second 
kind f22||, and p° = + p'^Y^^ the particle energy. 
Later on, the distribution parameter (3 will be identi- 
fied with the inverse thermodynamic (rest) temperature 
of the gas. The exact functional form of the spatial 
density g in Eq. pap is irrelevant, as long as g is nor- 
malizable (i.e., restricted to a finite spatial volume set 

V C M.^ in S). For simplicity, we may consider a spa- 
tially homogeneous gas enclosed in a stationary cubic box 

V = [—L/2, in E, corresponding to 



g{x) 



0, if cc ^ 



(3c) 



Here, V ~ is the S-simultaneously measured 

(Lebesgue) volume of V in S. 

The phase space PDF / is a Lorentz scalar Thus 
the current density and energy-momentum tensor 9^"^ 
can be constructed from / by: 



N 



9^^t,x) = N J fp^pr 



Concretely, for Eqs. ^ we have (j") = (p, 0) and 
I 0, A* / 



(4a) 
(4b) 

(5) 
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FIG. 1: The different hypersurfaces as defined in Eq. (|6]) 
and The worldlines of two particles are labeled by "PI" 

and "P2"; the worldlines of the container walls correspond 
to vertical lines at a; = —L/2 and x — L/2 (dashed lines), 
respectively. Assume a lab observer, resting at position a: = 
in E, takes a photograph of the system at the spacetime event 
S" with coordinates {t, x) — {^'', 0) in E. This photograph will 
reflect the state of the system along the backward-lightcone 
ff(<f). 



where = 3/?~^ -I- m Ki{Pm) / K2{Pm) is the mean 
energy per particle. One readily verifies that daj" = 0, 
whereas 9^6'^* = (3^^dig 7^ at the boundary of V, in 
agreement with Eq. ([T|). Confinement generates stress - 
the importance of this seemingly trivial statement shall 
be seen immediately. 

Isochronous state variables.- The traditional versions 
of relativistic thermodynamics 0, S, H, [13, [IH, El, [II 
are recovered from Eqs. (O-® by inserting isochronous 
spacetime hypersurfaces into Eq. To see this, con- 
sider an inertial frame S', moving at velocity w along 
the aj-'^-axis of the lab-frame E. An event S with coordi- 
nates ('C^,^) in S and (^"',^') in S' defines isochronous 
hypcrplanes 3(^°) and J'(C'°) in S and S', respectively, 
by 

3{e) := {{t,x)\t = e}, (6a) 
y^e) ■■= {{t',x')\t' = ^'^}. (6b) 

If E and S' are in relative motion, these hyperplanes dif- 
fer from each other, 3{£,'^) ^ 3'(^'°), see Fig.[T] Inserting 
= ^K*^] into Eq. we obtain the lab-isochronous 
energy- momentum vector 11 '^[3] in E: 



m3]) = N{{p°),0). 



(7) 



On the other hand, choosing S) = yields the E'- 

isochronous energy-momentum vector IC^p'] in S': 



'7((P°) + ^2/3-1), ^==0, 
0, A* > 1, 



(8a) 



where 7 := (1 — w^) 
mation with — w to Eq, 



Applying a Lorentz transfer- 
al), we find the corresponding 



3 



energy-momentum in E 

U^p'] = iV((p") ,-w/3-\0,0). (8b) 

Hence, the energy-momentum vectors ([7]) and (|8ap are 
not related by a Lorentz transformation. In fact, ICp] 
and U^p'] are connected by 

(u^'p']) = (u^p]) + iv/ri(o,w,o,o), (8c) 

reflecting the underlying hypersurface and observer ve- 
locity. As mentioned earlier, the difference between 
lt'^[3'] and 'U''[rJ'] arises because the energy- momentum 
tensor of a spatially confined gas is not conserved. It 
is also a reason for the existence of various temperature 
Lorentz transformation laws. 

Entropy.- Having at hand the state variables 'en- 
ergy' U*^ and 'momentum' U, one still needs 'entropy'. 
For a Jiittn er gas, one can define the entropy density 
four-current [2J, [2a| by (units fee = 1) 



^{t,x) 



-N 



d^p 



p'' fln{h\f). 



(9) 



The specific shape ([9]) is tightly linked to the exponential 
form of the Jiittner distribution ([3]). In fact, this com- 
bination (/, s'') is just one amongst several probabilistic 
models of thermodynamics; i.e., there exist other pair- 
ings, e.g., based on Renyi-type entropies, that yield con- 
sistent thermodynamic relations as well (26| . However, 
inserting ^ into we find 



^ ' ^ ^ [0, 1/ > 0. 



(10) 



Hence, the current PT?)) is stationary in S and satisfies 
the conservation law 



d^^s"" EE 0. 



(11) 



The associated thermodynamic entropy § is obtained by 
integrating over some space-like or light-like hyper- 
plane i3, yielding the Lorentz invariant quantity 



§[Sj] := [ da^ s^it.x). 



(12) 



Equation (fTT|) implies that the integral (fT2|) is the same 
for the hyperplanes 3(^°) and U'(^"^), 

§p] = §'p] = Sp'] = §'p']. (13) 

Thus, there is little or no room for controversy about 
the transformation laws of entropy in this example. The 
integral (fT2|) is most conveniently calculated along — 
3[£,^] in E, yielding 

§ = y" d^x s° = N ln{VZ/h^) + (iN (/) . (14) 

This can also be rewritten as 

§' = 7Vln(7y'Z/ft3) +/3U'°p]/7 (15) 
= N\n{-iV'Z/t?) + /37(U'"p'] + wU'^p']), 

where V' ~ V/^ denotes the Lorentz contracted (i.e., 
E'-simultaneously measured) volume [4l| . 



Einstein- Planck theory.- We are now ready to sum- 
marize the most common versions of relativistic thermo- 
dynamics. Planck Q and Einstein Q propose to use 
the E'-synchronous four- vector U'^p'] from Eq. ([5a|) as 
thermodynamic energy-momentum state variables. Fur- 
thermore, they choose to define heat Q'[J'] and, thus, 
temperature 7' in E' by the following postulated form of 
first law of thermodynamics [J^l 

dQ'p'] := T'dS' := dlt'°p'] - w'dU'^p'] + T'dl/', (16a) 

where the intensive variable w' = —w is the constant x'^- 
velocity of the gas (container) in E', and 3" the pressure. 
Considering the special case w' = first, we see that 
Eq. (|16ap is consistent with the second line of Eq. ([15]) 
upon identifying T = (i^^ and TV = Nf3~^; i.e., the 
parameter (3 of the Jiittner distribution equals the inverse 
rest temperature. Furthermore, for moving systems with 
w' 0, we find that thermodynamic quantities in E and 
E' are related by [H 

V' = y/7, 5" = T, §' = 8, (16b) 

lt'°p'] = 7(ll"p]+u;'2w) , (16c) 

It'^J'] ^ jw' (U°p] + yv) , (16d) 

T' = 7-iT=(l-ii;'2)i/2T, (i6e) 

i.e., within the Einstein-Planck formalism a moving body 
appears cooler [44|. Equations (|16p were criticized in a 
posthumously published paper by Ott [lOll and, later, also 
by van Kampen [§, and Landsberg [ill. 12 . 

Ott's and van Kampen's theory.- Ott 10[ and Van 



Kampen Q choose to formulate thermodynamic relations 
in the moving frame E' in terms of the E-isochronous 
energy-momentum vector U'^ [3] = A^j^U'^p]. They dif- 
fer, however, as to how heat and work should be de- 
fined. Van Kampen [23| replaces Planck's version of 
the first law, Eq. p6ap . by introducing a covariant ther- 
mal energy-momentum transfer four-vector Q'' via 



dQ^p] :=du^ p] -dyi^p]. 



(17) 



where, in the (lab) frame E, the non-thermal work vector 
yi^'p] is determined by (d^l'^p]) := {-?dV,0). Accord- 
ingly, in a moving frame E' one then finds dQ'^[3] = 
dlC' [3] — dyi'^ [3] , where by means of a Lorentz transfor- 
mation 

dll'^p] = w'" dlX"p], dyi'^ ^ -w'" 7dV. (18) 

Here, (w'^) ~ (7, 71/7', 0,0) denotes the velocity four- 
vector of the gas (container) in E'. While essentially 
agreeing on Eqs. pT[) . ([T8[) . and on the scalar character of 
entropy S' = §, van Kampen and Ott postulate different 
formulations of the second law, respectively. Specifically, 
Ott [l3| defines the temperature T' in E' via 



r d§' dQ'" = 7dQ" = 7 Ids, 



(19a) 



which implies the modified temperature transformation 
law El 



r = 7T=(l 



w 



^)-i/2T; 



(19b) 
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i.e., according to Ott's definition of heat and temperature 
a moving body appears hotter. Van Kampen [9| argues 
that the Eqs. p9|) are not well-suited if one wishes to 
describe heat and energy-momentum exchange between 
systems that move at different velocities (hetero-tachic 
processes) . To achieve a more convenient description, he 
proposes to characterize the heat transfer by means of a 
heat scalar Q' = Q, defined by % 2^ 

dQ' := -w'^AQ'" = -w^dQ" = dQ = dQ°. (20a) 

He then goes on to define temperature with respect to Q, 

T'dS' := dQ' = dQ = Td§, (20b) 

yielding yet another temperature transformation law: 

7' = T; (20c) 

i.e., according to van Kampen's definition a moving body 
appears neither hotter nor colder. Adopting this conven- 
tion, one can define an inverse temperature four-vector 
B' := T'^^ w' = T^^ w' and rewrite the second law in 

M H' P 

the compact covariant form 

d§' = -p'^dQ'f. (21) 

Discussion.- Evidently, whether a moving body ap- 
pears hotter or not depends solely on how one defines 
thermodynamic quantities. The formalisms of Ott [lol | 
and van Kampen are based on the same (lab-) 

isochronous hyperplane 3 - they merely differ in their 
respective temperature definitions p^ . By contrast, the 
Einstein-Planck theory 0, [Ml is based on an observer- 
dependent isochronous hyperplane 3'. While, in princi- 
ple, there is nothing wrong with this, a conceptual down- 
side of the latter approach is that the state variables en- 
ergy and momentum, when measured in different frames, 
are not connected by Lorentz transformations - or, put 
differently: To experimentally determine, e.g., the ener- 
gies 11° [3] and U"^ [3'] , two observers need to perform non- 
equivalent measurements [l8| . since measurements must 
be performed S-simultaneously in the first case, but E'- 
simultaneously in the second case. This might seem suf- 
ficient for regarding either Ott's or van Kampen's (more 
elegant) approach as preferable. However, before adopt- 
ing this point of view, it is worthwhile to ask: 

• Which formulation is feasible from an experimental 
point of view? 

• Which formalism provides a suitable conceptual 
framework for extensions to general relativity f2^ . 
\M]? 

Unfortunately, from an objective perspective, neither of 
the above proposals fulfils these criteria. The reason is 
that either formulation is based on simultaneously de- 
fined averages. On the one hand, this means that it is 
virtually impossible to directly measure the quantities 
appearing in the theory; e.g., in order to determine 11'^ [3] 



one would have to determine the velocities of the particles 
at time t = in S, which requires either superluminal 
information transport [Tsl or unrealistic efforts of trying 
to reconstruct isochronous velocity data from recorded 
trajectories. On the other hand, it is very difficult, if not 
impossible, to transfer the concept of global isochronicity 
to general relativity due to the absence of global incrtial 
frames in curved spacetimc. 

'Photographic' thermodynamics.- To overcome these 
drawbacks, we propose to define relativistic thermody- 
namic quantities by means of surface integrals along the 
backward-lightcone ^[S"], where S" is the event of the ob- 
servation, see Fig. [TJ This is motivated by the following 
facts: (i) A photograph taken by an observer at the 
event (o reflects the state of the system along the light- 
cone (ii) The hyperplane is a relativistically 
invariant object which is equally accessible for any in- 
crtial observer; i.e., if another observer 0', who moves 
relative to 0, takes a snapshot at the same event S", then 
the resulting picture will reflect the same state of the 
system - although the 'colors' will be different due to 
the Doppler effect caused by the observers' relative mo- 
tion [2^. (iii) The concept of lightcone-averaging can 
be easily extended to general relativity, (iv) In the non- 
relativistic limit c — > oo, the lightcone flattens so that 
photographic measurements become isochronous in any 
frame in this limit. Thus, lightcone averages appear to be 
the best-suited candidates if one wishes to characterize 
relativistic many-particle systems by means of nonlocally 
defined, macroscopic variables. 

Mathematically, the backward-lightcones in S 

and €.'[S'] in S' are given by 

e:(^) := { it,x)\t = e }, (22a) 
:= {{t',x')\t'^i''>-\x'-^'\}. (22b) 

Unlike the isochronous hyperplanes 3(^°) and 3'{£,"^), 
the lightcones describe the same set of spacetime events, 
= €'{£"). Fixing = £(<?) in Eq. Q, we flnd (see 
App.IH) 

1(°[G:] = iV(pO), (23a) 
= - J d'xj^^gix), (23b) 

Unlike U[3] and U[3'], the vector VL[€.] depends on the 
space coordinates ^ of the observation event S". Light- 
cones are Lorentz-invariant objects, implying that "[{''[C] 
and 'U'''[(£:'] are directly linked by a Lorentz transforma- 
tion, i.e., 'U''"[£'] = Ai^U''[€]. Moreover, for a spatially 
homogeneous Jiittner gas, it is straightforward to com- 
pute the entropy §[£] as [cf. Eq. (fT5|) ] 

§[£] Arin(UZ/ft3) +/3U°[£] 

= Nln{-fV'Z/h^) + /37(U'°[£] - w'U'^[€]) 

= S'M, (24) 

where, additionally, §[£] = §p] due to the conservation 
law (jlip . Thus, depending on which definition of heat we 
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choose, we again end up with cither Ott's or van Kam- 
pen's temperature transformation law. In our opinion, 
van Kampen's approach is more appealing as it defines 
temperature (similar to rest mass) as an intrinsic prop- 
erty of the thermodynamic system, whereas Ott's for- 
malism treats temperature as a dynamic quantity simi- 
lar to the zero-component of the energy-momentum four- 
vector. 

Observable consequences.- Since, unlike their 
isochronous counterparts, the state variables 1X''[£] are 
experimentally accessible, it is worthwhile to discuss 
implications for present and future astrophysical obser- 
vations. Let us assume that an idealized photograph, 
taken by an observer at S", encodes both the positions 
and velocities (e.g., from Doppler shifts) of a confined 
gas. If is at rest relative to the gas, then the mean 
values of the energy and momentum sampled from 
the photographic data will converge to given by 

Eq. Equation (|23ap implies that it does not matter 

for an observer at rest in S whether energy values are 
sampled from a photograph or from E-simultancously 
collected (i.e., reconstructed) data. 

The situation is different, when estimating the mean 
momentum from photographic data. Equation (j23bp im- 
plies that the lightcone-average depends on the observer 
position ^ [i^. A distinguished 'photographic center-of- 
mass' position in S can defined by 

^'M||=l = ' z = 1,2,3. (25) 

For example, if g is symmetric with respect to the origin 
of E, then = 0. This would correspond to a lightconc 
as drawn in Fig. [1] In this (and only this) case, we find 
= i(j"U'°[£] and. thus, lightconc thermodynamics 
reduces to the Ott-van Kampen formalism. 

To illustrate how 1X'[£] generally depends on the ob- 
server's position, let us consider a gas with density pro- 
file ([3c]). For a stationary observer at a position ^ far 
away from V, we can approximate |a; — ^| ~ |^| in 
Eq. ((23b)) . yielding ^ 

U'[t\ = - ili^T-. (26) 

Hence, a distant observer who naively estimates 'U*[£] 
from photographic data could erroneously conclude that 
the gas is moving away with a momentum vector propor- 
tional to the temperature. Reinstating constants c and 
fee, this relativistic effect becomes neglible if mc^ ^ ks'^, 
but - given the rapid improvement of telescopes and 
spectrographs [30| - it should be taken into considera- 
tion when estimating the velocities of astrophysical ob- 
jects from photographs in the future. In particular, since 
Lorentz transformations mix energy and spatial momen- 
tum components, for a moving observer 0' both 'U'°[£] 
and it'^fC] will be affected, see Fig. |4l In principle, simi- 
lar phenomena arise whenever one is limited to photo- 
graphic observations of partial components in distant 
compound systems (e.g., the gas in a galaxy), if the 




FIG. 2: The distribution of the particle energy p in the rest 
frame of the gas for various values of the mean particle energy 
(p"): The distribution measured in our simulations (symbols) 
is compared to the Jiittner distribution ()3b[l (lines). 



energy-momentum tensor of this partial component is not 
conserved. At present, it is an open problem whether or 
not these effects may even be amplified in curved space- 
time geometries. 

Numerical results. - The preceding theoretical consid- 
erations can be illustrated by (l-|-3)-dimcnsional rela- 
tivistic many-body simulations. Compared with the non- 
relativistic case, simulations of relativistic many-particle 
systems are more difficult because particle collisions can- 
not be modeled by simple interaction potentials any- 
more dll, [13 . Generalizin g re cently proposed lower- 
dimensional algorithms pll . l33l | , our computer experi- 
ments are based on hard-sphere-type interactions in the 
two-particle center-of-mass frame (see App. [Cl details). 
This model is fully relativistic in the low-to-intermediate 
density regime [4^. Figures [JHl] depict results of simula- 
tions with N — 1000 particles. Initially, all particles are 
randomly distributed in a cubic box with same energy p", 
but random velocity directions. After a few collisions per 
particle, the energy distribution relaxes to the Jiittner 
distribution (j3bp , see Fig. [21 which confirms that our col- 
lision algorithm works correctly in this density regime. 
The thermodynamic energy-momentum vector 'U''^[i3] is 
determined by recording each particles' momentum as its 
trajectory passes through the corresponding hyperplane. 
We first consider an observer who is at rest relative 
to the gas. As predicted by Eq. ([26]) we find that, if the 
location of deviates from the center of the box, a photo 
made by yields a non-zero momentum IX' cx see 
blue line in Fig. [3| For comparison. Fig. [4| shows the 
results for a moving observer 0' (speed w), obtained by 
isochronous sampling along different hyperplanes y{w), 
or photographic sampling along the backward-lightcones 
€((?). Again, as predicted by the theory, the resulting 
overall energy-momentum does not only depend on the 
observer velocity, but also on the underlying hypersurface 
and, in particular, on the observer's position. Although 
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FIG. 3: The lightcone momentum U^[£] as a function of gas 
temperature for an observer resting either at the center of 
the gas container, ^ = 0, or far away from it, ^ — > (— oo, 0, 0). 
Symbols indicate values measured in the simulations, and 
lines the theoretical predicted curves. 



our study still neglects quantum processes and gravity, 
which play an important role in real astrophysical sys- 
tems, the results suggest that one needs to be very care- 
ful when reconstructing the velocities of hot, relativistic 
objects from photographic measurements. 

Conclusions.- Sometimes, discussions of relativistic 
thermodynamics start by postulating a set of macro- 
scopic state variables, whose thermodynamics relations 
(and Lorentz transformations laws) are subsequently de- 
duced by plausibility considerations. Unfortunately, this 
approach - although quite successful in nonrelativistic 
physics - is intrinsically limited in a relativistic frame- 
work, as it conceals the actual source for conceptual diffi- 
culties, namely, the nonlocal character of thermodynamic 
quantities. The above analysis may provide guidance for 
constructing consistent relativistic thermodynamic the- 
ories for more complex systems, e.g., based on (non-)- 
conserved tensor densities as derivable from classical or 
quantum Lagrangians. Care is required when integrat- 
ing such densities to obtain global thermodynamic state 
variables, since conservation laws may be violated due 
to confinement, so that averages may vary depending on 
the underlying hyperplane(s). Within a conceptually sat- 
isfying and experimentally feasible framework, thermo- 
statistical averaging procedures should be defined over 
lightconcs rather than isochronous hypersurfaces. To 
put it somewhat provocatively, the isochronous defini- 
tion of nonlocal quantities, as adopted in traditional for- 
mulations of relativistic thermodynamics, can be viewed 
as a relic of our accustomed nonrelativistic 'thinking'. 
With regard to present and future astrophysical obser- 
vations, it will be important to better understand how 
the temperature-dependent, apparent drift effects due 
to lightcone-averaging (i.e., photographic measurements) 
become modified in curved spacetime, because this might 
affect velocity estimates for astronomical objects which 
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FIG. 4: Observed energy U'" (a) and momentum U'^ (b) of 
the gas as a function of the observer speed w along the x^- 
axis, either sampled from an isochronous hyperplane 3', or a 
lightcone The mean particle energy is (p") = 3m. The 
observer using a lightcone is either at the center of the gas 
container (^ — 0), far behind the container (^^ — » — oo), or 
ahead of the container (^^ -fcjo). Simulation results are 
indicated by symbols, the theoretical prediction by lines. 



are pivotal for our understanding of the cosmological evo- 
lution M- 



APPENDIX A: NOTATION 



We adopt units such that speed of light c = 1 and 
Boltzmann constant A:b = 1, and the metric convention 
iVfiiy) = diag(— 1, 1,1,1) = {v'fj.iy)- Einstein's sum rule is 
applied throughout. If an event has coordinates {^^) = 
(C",f ) = (?°,^^'?^C^) = (t,^) in the inertial spacetime 
frame E, then its coordinates in another frame S', moving 
at constant relative velocity w along the x^-axis of S, are 
given by (e"^) - {i{e ~ <'),7(e' - <°),C',^') with 
7 



{1 - u)2)-i/2. In short, = ^'"uS,", where (A'^^) is 



the corresponding Lorentz transformation matrix. 
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APPENDIX B: SURFACE INTEGRALS IN 
MINKOWSKI SPACETIME 



APPENDIX C: NUMERICAL SIMULATIONS 



We wish to calculate the integral 



J S^ 



(Bl) 



where is a three-dimensional hyperplane in the (1 + 3)- 
dimensional Minkowski frame E. If O^"^ " is a tensor 
of rank n then S"'^" [i3] has rank (n — 1). Considering 
Cartesian spacetime coordinates, the surface element du^ 
may be expressed in terms of the alternating differential 
form [11 



dff^ = (3!)"^e^„;3^d.T" A da;*^ A da;'', 



(B2) 



where e^a/ij is the Levi-Cevita tensor |35| . and da;" A 
Ax^ = — da;'^ A da;" the antisymmetric product. With re- 
gard to thermodynamics, we are interested in integrating 
over space-like or time-like surfaces Sj given in the form 
= t = g{x). Examples are the isochronous hyperplane 
from Eq. dH) and the hghtcone from Eq. ((22)) . 
Given the function we may write da;*' = dig da;*. In- 
serting this expression into Eq. (jBip yields 



S" - = / dx^Adx^Ada;^ [6i°"- - 6i*" -] 

d?x [0°"-(.g(a;),a;) - {d,g) O^'^- {g{x),x)\ 



In particular, for the isochronous hyperplane from 
Eq. ([HI), we have g{x) = and dig = in S, leading to 



(B3) 



For the lightcone with dig — — (x* — ^*)/|a; — we 
find 

T-\^\ - yd^x + 

f^0"-(^°-|a;-^|,a.)l. (B4) 



Fully relativistic iV-body simulations would require to 
also compute the interaction fields generated by parti- 
cles, which is numerically expensive. For dilute gases 
with short-range interactions one can obtain reliable re- 
sults by considering simplified models based on quasi- 
elastic collisions. In our computer experiments, we simu- 
late a three-dimensional gas of relativistic hard spheres in 
a cubic box. In a particle-particle collision, momentum 
is transferred instantaneously at the moment of closest 
encounter by taking into account the relativistic energy- 
momentum conservation laws. 

The tasks during a simulation time-step are 
(i) Determine the times/distances of all particle pairs at 
their closest encounter, (ii) Advance all particles to the 
next collision time, (iii) Transfer momentum between 
the colliding particles, (iv) Record particle energies and 
momenta when the particles are reflected at the walls 
(e.g., to measure the pressure on the boundaries) or their 
trajectories pass an observer hypersurface. 



Our simulations show that details of the momentum 
transfer mechanism (e.g., the differential cross-sections) 
do not affect the stationary momentum distribution. It 
is, however, important to employ a Lorentz- invariant col- 
lision criterion (we use the minimum distance of particles 
in the two-particle centcr-of-mass frame). Non- invariant 
criteria (e.g., the minimum intcr-particlc distance in the 
lab frame) may lead to deviations from the Jiittner dis- 
tribution. 



The most time-consuming task is to determine the 
closest-encounter times/distances for all particle pairs. A 
huge speed-up can be achieved by considering only close 
particle pairs, using a hash table based on a partition of 
the simulation box into subcubes. With this method one 
can efficiently simulate 10'^ particles and 10® collisions in 
30 min on a desktop pc. 
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More precisely, one should write V' = V^'[3'] and V — 
V[3] in order to reflect how volume is measured (defined) 
in either frame. 
Cf. Eq. (23) in [g. 

The scalar transformation law of the pressure is consis- 
tent with the Lorentz transformation laws of force and 
area [^. 

It seems that, in the later stages of his life, Einstein 

changed his opinion about the correct transformation 

laws of thermodynamic quantities, favoring formulas 

which were later independently derived by Ott [l3| and 

Arzelies [H, [131 ; see Liu [H, [13 

See also Eddington [i^l and Arzelies [36l |. 

Averages in the lab frame do not depend on the specific 

value of the time coordinate if the PDF is stationary 

in this frame. 

Equations (|23|l and (|26|) are consistent with Eq. ((8)1, as 
can be seen by letting w ^ 1 in H8b[l and 3> 0, = 
= in (|26|l. 

Generally, hard-sphere models become 'nonrelativistic' at 
very high densities (liquid/solid regime). 



